Wave turbulence on the surface of a ferrofluid submitted to a magnetic field 
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We report the observation of wave turbulence on the surface of a ferrofluid mechanically forced and 
submitted to a static normal magnetic field. We show that magnetic surface waves arise only above 
a critical field. The power spectrum of their amplitudes displays a frequency-power law leading to 
the observation of a magnetic wave turbulence regime which is experimentally shown to involve a 
4- wave interaction process. The existence of the regimes of gravity, magnetic and capillary wave 
turbulence is reported in the phase space parameters as well as a triple point of coexistence of these 
three regimes. Most of these features are understood using dimensional analysis or the dispersion 
relation of the ferrohydrodynamics surface waves. 

PACS numbers: 47.35.Tv,47.65.Cb,47.27.-i 



Wave turbulence is an out-of-equilibrium state where 
waves interact with each other nonlinearly through TV- 
wave resonance process. The archetype of wave turbu- 
lence is the random state of ocean surface waves, but it 
appears in various systems: capillary waves plasma 
waves in solar winds, atmospheric waves, optical waves, 
and elastic waves on thin plates 0]. Recent laboratory 
experiments of wave turbulence have shown new obser- 
vations such as intermittency j^, fluctuations of the en- 
ergy flux 5] , and finite size effect of the system 0, Q . 
Some of these phenomenon have recently been considered 
theoretically [7|. Wave turbulence theory allows to ana- 
lytically derive stationary solutions for the wave energy 
spectrum as a power law of frequency or wave number [3| . 
The spectrum exponent and the number N of resonant 
waves depend on both the wave dispersion relation and 
the dominant nonlinear interaction. Several theoretical 
questions are open, notably about the validity domain 
of the theory [Sl], and the possible existence of solutions 
for non dispersive systems Q. In this context, finding 
an experimental system where the dispersion relation of 
the waves could be tuned by the operator should be of 
primary interest to test the wave turbulence theory. 

A ferrofluid is a suspension of nanomctric ferromag- 
netic particles diluted in a liquid displaying striking prop- 
erties: the Rosensweig instability 10| , the labyrinthine 
instability, magnetic levitation [ll'|. In contrast with 
usual liquids, the dispersion relation of surface waves on 
a ferrofluid displays a minimum which depends on the 
amplitude of the applied magnetic field 12,[l3|. Thus, 
one can easily tune the dispersion relation of surface wave 
from a dispersive to a non dispersive one with just one 
single control parameter. To our knowledge, no exper- 
imental observation of wave turbulence on a magnetic 
fluid has been reported. Here, we study the wave turbu- 
lence on the surface of a ferrofluid submitted to a normal 
magnetic field. We observe for the first time a regime of 
magnetic wave turbulence. We characterize this regime 
by measuring the power spectrum and distribution of the 



FIG. 1: Experimental setup 



magnetic wave amplitude. 

The experimental setup is shown in Fig. [TJ It consists 
of a cylindrical container, 12 cm in inner diameter and 4 
cm in depth, filled with a ferrofluid up to a depth h — 2 
cm. The ferrofluid used is a ionic aqueous suspension 
synthesized with 8.5% by volume of maghemite particles 
(Fe203 ; 7.6 ±0.36 nm in diameter) [14]. The properties of 
this magnetic fluid are: density, p = 1324 kg/m'^, surface 
tension, 7 = 59 x 10^^ N/m, initial magnetic suscepti- 
bility, Xi — 0.69, magnetic saturation Mgat — 16.9 x 10'^ 
A/m, and estimated dynamic viscosity 1.2 x 10"'^ Ns/m^. 
The container is placed between two horizontal coaxial 
coils, 25 cm (resp. 50 cm) in inner (resp. in outer) diam- 
eter, 7 cm far apart. A DC current is supplied to the coils 
in series by a power supply (50 V/35 A), the coils being 
cooled with water circulation. The vertical magnetic in- 
duction generated is 99% homogeneous in the horizontal 
plane [I4I, and is up to 780 G. It is measured by a Hall 
probe located in the center near the surface of the con- 
tainer. Surface waves are generated on the ferrofluid by 
the horizontal motion of a rectangular plunging Teflon 
wave maker driven by an electromagnetic vibration ex- 
citer. The wave maker is driven with low-frequency ran- 
dom vibrations (typically from 1 to 5 Hz). The amplitude 
of the surface waves ri{t) at a given location is measured 
by a capacitive wire gauge (plunging perpendicularly to 
the fluid at rest) with a 7.1 mm/V sensitivity [3]. ri{t) is 
recorded by means of an acquisition card with a 4 kHz 
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sampling rate, low-pass filtered at 1 KHz during 300 s, 
leading to 1.2 x 10^ points recorded. 

In the deep fluid approximation, the dispersion rela- 
tion of linear inviscid surface waves on a magnetic fluid 
submitted to a magnetic induction B perpendicular to 
its surface, reads 
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where lo is the wave pulsation, k its wavenumber, g = 
9.81 m/s^ the acceleration of the gravity, /io = 47r x 
10^^ H/m the magnetic permeability of the vacuum, and 
/[x] = x^/[(2+x)(l+x)]- X is the magnetic susceptibility 
of the ferrofluid which depends on the applied magnetic 
field, iJ, through Langevin's classical theory [3] 
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where C{x) = coth (a;) — 1/x, and thus on the magnetic 
induction, B, through an implicit equation since 



B^^loil + x)H 
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For B = 0, the dispersion relation of Eq. ([T]) is 
monotonous, and is dominated by the gravity waves at 
small k, and by the capillary waves at large k. When B 
is increased, the quadratic term —B^k^ increases, and 
the dispersion relation becomes non monotonous: an 
inflection point appears at _B = 0.9'iBc, then a mini- 
mum which leads, at i? = B^, to the Rosensweig in- 
stability [lH. This stationary instability (a hexagonal 
pattern of peaks on the ferrofluid surface) occurs when 
a'^(fc) becomes negative, that is for a critical induction 
-B^ = 2iJ,n J pg^/ f\\( He)] where xi^c) is determined us- 
ing Eqs. ^ and (0]) [XO]. This leads to the theoreti- 
cal value of the critical induction Be = 292.3 G for the 
Rosensweig instability of our ferrofluid. When B is slowly 
increased, a direct visualization of the surface leads to an 
experimental value of Be = 294 ± 2 G which is close to 
the above expected value. 

The power spectrum of the wave amplitude on the sur- 
face of the ferrofluid is shown in Fig. [2] for different ap- 
plied magnetic induction B. For B = (see the inset of 
Fig. [2]) , it displays similar results than those found with a 
usual fluid : two power laws corresponding to the grav- 
ity and capillary turbulence wave regimes. The capillary 
regime is found to scale as J~2.9±o.i jj-^ good agreement 
with the prediction of weak turbulence theory in /"^^^^ 
[l6| . and the gravity regime is found in f~^'^. The expo- 
nent of the gravity cascade being known to depend on the 
forcing parameter 0,0] (in contrast with the theoretical 
prediction ~ 1121 )j it is thus only fitted to measure 
the crossover frequency fgc between gravity and capillary 
regimes. As previously reported with usual fluids [2], fgc 
is also found here to decrease (from 26.6, 21 to 17.2 ±0.3 
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FIG. 2: Power spectrum of rj{t) for two values of B. In- 
set: 5 = 0: Gravity and capillary wave turbulence regimes. 
Dashed lines have slopes -4.6 and -2.9. Crossover: fgc — 20 
Hz. Main: B — O.QBc'. Magneto-capillary wave turbulence. 
Dashed line has slope -3.3. Crossover: fgc — 5 Hz. Forcing 
parameters: 1-5 Hz. Be — 294 G. 



Hz for a random forcing of 1 - 6 Hz, 1-5 Hz, and 1-4 Hz). 
The expected value is given by fgc — ^y^2g/lc — 15.2 Hz 
where Ic = \/l/{pg) 0]- Consequently, in order to study 
the crossover frequency dependence with the magnetic 
induction B, one has to rescale it by its value at i? = 
such as fiB) = fgc{0) X f{B)/f{Q) (see below). 

For _B 7^ 0, the power spectrum of wave amplitudes 
shows two striking new results. As shown in Fig. [51 the 
crossover frequency is strongly decreased down to roughly 
5 Hz, and a power law in f^^'^ appears in roughly all 
the accessible frequency range. Let us try first to under- 
stand this latter observation. The power spectrum oi'q{t) 
can be derived by dimensional analysis for the gravit 
and capillary wave turbulence regimes, respectively as 
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where e is the energy flux [dimension (L/T)^]. Since grav- 
ity wave turbulence is a 4-wave interaction process, and 
capillary waves a 3- wave one, it gives the dependence of 
the energy flux exponent in for a iV-wave pro- 

cess One can also derive dimensionally the power 
spectrum for the magnetic wave turbulence regime as 
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In contrast with the above dispersive systems, this —3 
frequency exponent does not depend on the energy flux 
exponent a, that is on the number N of resonant waves. 
The frequency exponent predictions for the magnetic and 
capillary regimes, respectively —3 and —17/6 ~ —2.8 can 
not be distinguished experimentally within our experi- 
mental accuracy. This could explain that only one single 
slope is observed on Fig. [2] which thus corresponds to a 
"magneto-capillary" wave turbulence regime. 
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FIG. 3: Exponent of the magneto-capillary spectrum as a 
function of the dimensionless magnetic induction B/Bc- In- 
set: Amplitude of the power spectrum (averaged between 
8 and 18 Hz) as a function of B/Bc- Forcing parameter: 
1 < / < 4 Hz. Be = 294 G. 



FIG. 4: Rescaled crossover frequencies, f{B) = /gc(0) x 
f{B)/f{0) as a function of B for different frequency band- 
widths of the random forcing: (x) 1 to 4 Hz, (o) 1 to 5 Hz, 
and (+ or o) 1 to 6 Hz. Theoretical curves fgc, fgm and fmc 
are respectively from Eqs. ((5}, (|6]) and ((7)|. The triple point 
{ft = 10.8 Hz , Bt/Bc = 0.65) is from Eq. JH}. 



Figure[3]shows the frequency exponent of the magneto- 
capiUary spectrum when B is increased. For small B, 
capillary waves are dominant and the exponent is found 
roughly constant ~ —2.9 in good agreement with the cap- 
illary prediction. When B/Bc > 0.65, magnetic waves 
becomes dominant (see below), and the exponent begins 
to slightly increase with B up to -3.1 in rough agreement 
with the -3 prediction of Eq. (g]). For B / B^ > 1, the 
Rosensweig instability occurs, and the spectrum expo- 
nent strongly changes with B. This could be attributed 
to the growing of the hexagonal pattern. The inset of 
Fig. [3] shows the power spectrum amplitude (averaged 
between 8 and 18 Hz) as a function of B. The spectrum 
amplitude is found roughly constant when B is increased 
till B/Bc — 0.65, onset of the magnetic waves. Above 
this value, the spectrum amplitude increases roughly lin- 
early with B up to B/Bc = 1 where the instability oc- 
curs. Using Eq. (|4]), one can thus deduce from this linear 
dependence that 2 — 3q; = 1, and thus a = 1/3. This 
experimentally shows that the magnetic wave turbulence 
observed involves a 4-wave interaction process. 

The crossover frequency between the gravity and 
magneto-capillary regimes decreases when B is increased 
(see Figure [5]). Figure [3] then shows the evolution of 
the rescaled crossover frequency f{B) (see above) as a 
function of B for 3 different frequency bandwidths of the 
random forcing. When B is increased, / is found to de- 
crease with the same law whatever the forcing frequency. 
Beyond B = Be, f is roughly of the same order than the 
upper frequency of the forcing, and consequently can not 
be measured anymore. However, when B > Be and for 
the 1-6 Hz forcing, the power spectrum displays three 
power laws (not shown here) showing a magneto-capillary 
crossover. This new crossover frequency is reported in 
Fig. m with (o)-symbols without rescaling. For lower fre- 



quency bandwidths of vibration, the slope breaking is too 
small to allow an accurate measurement of the transition 
between the magnetic and capillary regime. 

The evolutions of these crossover frequencies with B 
are described as follows. Whatever B, Eq. (jT]) is dom- 
inated, at small k, by the linear term (gravity waves) 
and, at high k, by the cubic term (capillary waves). One 
can assume that the quadratic term (magnetic waves) 
dominates when it is greater than the linear and cubic 
terms. This arises when f[x]B'^ > ^^o^/pgl, that is, us- 
ing Eqs. ^ and ^ and the ferrofluid properties, when 
B > 0.65i?c- Thus, when B < O.Q5Bc, no wavelength 
exists for which the magnetic term dominates in Eq. ([T]). 
When B > 0.65i3c, magnetic waves exist on a range of 
wavelength between the gravity and capillary ones. This 
explains why both the spectrum exponent and its am- 
plitude shown in Fig. [3] change for B ~ 0.65i3c- Sur- 
prisingly, this critical magnetic induction has never been 
reported previously. The crossover frequencies between 
the gravity, the magnetic and the capillary regimes are 
derived by balancing the dispersion relation terms each 
to each. For the gravity-capillary transition, one balances 
the first and the third terms of the second-hand of Eq. 
([1]), that is gkgc = {'j/p)kgc, thus for kgc — \/ pg/^ which 
substituted into Eq. ([T]) gives 

^sc = 2W , for f[x]B < poy/pgj . (5) 

V '^0'^ 

Similarly, by balancing the first and second terms, the 
gravity-magnetic crossover frequency reads 
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Finally, by balancing the second and the third terms, the 
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FIG. 5: Probability density functions of wave amplitude, ri, 
for different values of the dimensionless magnetic induction, 
from B/Bc = 0, 0.3, 0.54, 0.77, 0.92 to 1.2 (see the arrows). 
Forcing parameter: 1 < / < 4 Hz. Bottom inset: Same 
PDFs displayed using the reduced variable rj/ar^. Gaussian 
fit (dashed line). Top inset: Standard deviation of wave am- 
plitude an as a function of B/Bc 



magneto-capillary crossover frequency reads 
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Using Eqs. ^ and ((3]) and the ferrofluid properties, these 
crossover frequency curves are plotted in Fig. [4] as a func- 
tion of B/Bc and show a good agreement with the exper- 
imental data. This plot also shows the range of existence 
of magnetic waves for B/Bc > 0.65 and of a triple point 
corresponding to the coexistence of the three domains 
(gravity, magnetic and capillary ones) . It can be derived 
by balancing the three terms of Eq. ([T]) , which leads to 
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^nd B^ = ^^oVP9l/nx(Ht)] , (8) 



corresponding to ft = 10.8 Hz and Bt/Bc = 0.65 in good 
agreement with the data of Fig. |4l 

Finally, the probability density functions (PDFs) of 
the wave amplitudes are shown in Fig. [5] for different 
value of B/Bc, at high enough amplitude of forcing. For 
5 = 0, the PDF is asymmetrical due to the strong steep- 
ness of the waves as for usual fluids [2|. This means 
that the deep through are rare, whereas high crests are 
much more probable, thus showing the nonlinear nature 
of the wave interactions in a wave turbulence regime. 
The PDF asymmetry is enhanced when B is increased. 
Note also that the most probable value of the PDF is 
more and more negative as B increases, although its 
mean value {if) is zero as expected. The bottom inset 
of Fig. [5] shows that all these PDFs normalized to its 
standard deviation value, u^, roughly collapse on a single 



non Gaussian distribution. This means that the distri- 
bution depends only on . The top inset of Fig. [5] then 
shows the evolution of function of B/Bc- For 

B/Bc < 0.63, the rms value amplitude of the waves does 
not depend on the magnetic induction. Note that this 
value is very close to the above predicted onset of mag- 
netic waves Bt/Bc = 0.65. When this onset is exceed, 
increases roughly linearly with B up to the occurrence of 
the Rosensweig instability at Be- 
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